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We discuss the application of quantum-mechanical supersymmetry to 
particle traps. The supersymmetric-partner wave functions may be 
used to describe a valence fermion in a trap system with an isotropic 
harmonic-oscillator potential. Interactions with the core are incorpo- 
rated analytically. The close similarity of this approach to the applica- 
tion of supersymmetry in atomic systems is made explicit by means of 
a radial mapping between the two systems. 
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Supersymmetry has been an active research area for well over two decades. De- 
spite this, very few physical supersymmetries are experimentally known. One is the 



appearance of an effective radial potential in the context of atomic systems [13], [15 
We discuss a possible further application in the context of particle traps. More details 
are given in the two references listed under our names. 

The isotropic harmonic oscillator has a radial equation admitting a supersymmet- 
ric partner. The physical implications of this mathematical fact may be investigated 
using trap systems for which an isotropic potential can be established. In the Ioffe- 



Pritchard trap [jTT], [12[ and the time-averaged orbiting-potential (TOP) trap |23 ], this 
condition may be satisfied. Both are neutral-particle traps that use the interaction 
of the magnetic dipole moment ft of the particle with a confining magnetic field. The 
traps select all dipoles aligned opposite to the direction of the magnetic field and 
draw them into the region of weakest field at the center. The Ioffe-Pritchard trap is 
purely magnetostatic, with the field provided by two coils and four linear conductors. 
The TOP trap comprises six coils, four of which have alternating currents, creating 
a high-frequency rotating magnetic field. An averaging procedure removes the time 
dependence, yielding an effective magnetostatic potential. With suitable choices of 
currents in the conductors of these systems |I|, |T9[ , it is possible to ensure isotropy of 
the potential energy near the center of the trap, 

U(r) = /2B Q (l + r 2 /r { 

Here, r is a length characteristic of the trap system, and B ^ is the magnitude 
of the magnetic field at the center of the trap. The radial wave functions describing 
a single trapped dipole may be expressed in terms of the generalized Laguerre poly- 
nomials Ltf (z). Ignoring a factor of r that removes the first-order derivative in the 
differential equation, the functions are 

W N , L (r) = C N , L (r/r ) L+1 exp (-r*/2ri)L%£% (r 2 /rl) . (2) 

Here, L = 0, 1, 2, ... is the angular momentum, and N = L, L + 2, L + 4, . . . is the 
principal quantum number. Normalization is ensured via the constant Cn,l- The full 



solutions of the harmonic oscillator \N, L, M) behave as 

W N , L (r)Y LtM (9,<!>) , (3) 

where the azimuthal quantum number M takes the usual values. 



Details of supersymmetric quantum mechanics may be found elsewhere [22, 26 



13, 1C[. We use the term bosonic sector to refer to the given radial system with fixed 
angular momentum L and spectrum shifted to have zero lowest-state energy. The 
term fermionic sector is used for the partner system. 

One of our objectives is to regard the radial-equation fermionic sector as providing 
an effective potential for an excited valence particle in a trap. It experiences not only 
the trapping potential, but also interactions with a core of other particles. Before 
motivating this application of supersymmetry further, let us consider the case of a 
valence particle with angular momentum L = that is excluded from occupying 
levels below N = 2 by a filled core. This exclusion can only hold for fermions, to 
which we restrict ourselves here. If the interactions between the trapped dipoles are 
small compared with the natural spacing of the energies in the oscillator, the number 
of particles in the core may be found by counting the levels in the single-particle 
bosonic system lying below that of the valence fermion. There are four in this case. 
One is the ground state, \N = 0, L = 0, M = 0), and the other three can be labelled 
as \N = 1, L = 1, M = 0, ±1). Similarly, the core would have 20 fermions if the L = 
valence fermion was restricted to N > 4. 

The physical interpretation of the fermionic sector is motivated by several ob- 
servations. It is well known that the fermionic sector is degenerate with the bosonic 
sector except for the lowest bosonic state, which has no corresponding fermionic state. 
One may visualize a situation in which physically such a corresponding state exists, 
but is inaccessible to the valence fermion because it is occupied by a core fermion. 
The absence of a zero-energy state in the fermionic sector of the mathematical for- 
malism reflects this physical picture. The core need not be occupied by only one 
fermion, since there might be others of different angular momenta. Adopting this 
interpretation, the Pauli principle is seen to underlie radial supersymmetry, and it 
becomes natural to interpret the fermionic sector as describing a valence fermion in 



a multifermion system. For L = 0, the effective radial potential for the fermionic 
sector differs from the corresponding one for the bosonic sector by an expression that 
includes the term h 2 /mr 2 where m is the valence fermion mass. This additional 
repulsion is thus consistent with the meaning of the Pauli principle. 

So, the effect of the supersymmetry is to fill an inner core with fermions. This 
procedure must leave the angular momentum of the valence particle unchanged. We 
therefore construct the full three-dimensional wave functions from the product of the 
fermionic radial wave functions Wjv 8 -i,.&+i( r ) an d the same spherical harmonics as in 
the bosonic sector (§): 

W^_ 1 , L+1 (r)Y LiM (^,0) , (4) 

where N s = L + 2, L + 4, L + 6, ... is the principal quantum number. Even though the 
constant L in the radial function appears to have been shifted, the angular momentum 
of the system is defined by the spherical harmonics and is unchanged. These full wave 
functions (f|) for the fermionic sector differ from the full wave functions (|3]) for the 
three-dimensional isotropic harmonic oscillator. 

An alternative description for multifermion traps could account for the filled core 
by using the standard harmonic-oscillator solutions (0) but requiring the values of the 
principal quantum number for the valence fermion to exclude numbers corresponding 
to the filled core. However, the valence fermion would then be described by an 
incomplete set of states and the lowest valence state would have too many nodes. Both 
of these drawbacks are absent for the fermionic functions obtained via supersymmetry, 
which form a complete orthonormalizable set of states and for which the lowest state 
has no nodes. They therefore resemble solutions for other conventional bound systems 
in quantum mechanics. 

The number of fermions in the core of a particular trap described by the functions 
(f|) depends on the angular momentum of the valence fermion. For L = 0, the princi- 
pal quantum number takes values N s — 2, 4, 6, . . . and by the reasoning considered 
above, this trap has four core fermions. The filling of another shell in the core can 
be accomplished by shifting the fermionic-sector spectrum to have zero lowest-energy 
state, thereby treating it as a new bosonic sector. A new fermionic sector is then 



obtained via the usual supersymmetry procedure. This describes a trap with 20 core 
fermions and a valence fermion with zero angular momentum. Further iterations of 
this procedure fills further shells, giving cores with 56, 120, 220, . . . fermions. For va- 
lence fermions with L — 1, the concept is the same, and there are 1, 10, 35, 84, . . . core 
particles. General formulae for these sequences can be obtained |T!J. They assume 
only one spin orientation since the dipoles in the Ioffe-Pritchard and TOP traps are 
oriented against the magnetic field. The formulae differ for other trap systems. 

Although it accounts for the Pauli principle, radial supersymmetry ignores inter- 
actions between the valence fermion and the core. We discuss one method of in- 
corporating interactions developed in analogy with analytical supersymmetry-based 
quantum-defect theory for atomic systems fi~6| . The modifications in the oscillator 
energy spectrum due to interactions plausibly generate the form 

E N * = fj,B + huj (N* + 3/2) , (5) 

for ujq = (2[xBq/ mr 2 ,) 1 ' 2 . The modified eigenvalues are thus incorporated via a shifted 
principal quantum number N*. Denoting the shift by A = A(N,L) and including 
also an integral shift / = I(L), we define N* = N + I — A. Equivalently, with 
N s = N + 21, we write N* = N s — I — A. If in addition we shift the angular 
momentum, L* = L + I — A, and add the effective potential 



Ti 2 L*(L* + 1)-L(L + 1) 
2m r 2 



VW(r) = ;Sr 1 2 V ; + M# - **) (6) 



to the differential operator in the Schrodinger equation for (|2]), analytical radial wave 
functions that correspond to the modified eigenspectrum (Kj) are obtained. They 
may be expressed in terms of the functional form (0) as Wjv* ) £,*(r). This analytical 
defect theory extends the radial supersymmetry quite naturally, and if the defects are 
switched off appropriately the exact bosonic and noninteracting fermionic sectors are 
recovered. 

The application of supersymmetry in traps closely follows the application of su- 
persymmetry in multi-electron atoms and in ions Jl4|, |j~5| . The similarity in the two 
applications is more than a mathematical parallel and can be made explicit in the 
form of a mapping. The existence of a natural correspondence between the radial 



three-dimensional Coulomb problem and the radial harmonic oscillator in two or four 



dimensions was first noted by Schrodinger more than fifty years ago p5[, and since 



then has received much attention [gj, g, [u], [Lq] . We consider Coulomb dimensions 
d > 1, to avoid normalization issues associated with the one-dimensional case, and 
oscillator dimensions D > 1. For these arbitrary-dimensional cases, the radial equa- 



tions may still be separated [^TJ, and we write the radial solutions as Wc^^{r) for the 
oscillator and Wd t n,i{ r ) fc> r the Coulomb case. Details of these functions may be found 
elsewhere [0. As an example, Wo= , i,L, v( r ) is identical to @). We adopt lower-case 



symbols for the Coulomb system and upper-case symbols for the oscillator systems, 
with an exception made for the oscillator radial variable r. In the Coulomb system, 
this convention gives angular momentum I and principal quantum number n. 

A natural mapping between these two radial systems exists subject to certain con- 
ditions on the dimensions, the angular momenta, and the principal quantum numbers. 
The relationship between the radial wave functions is [17] 



W D ^ L (r) = K d ^ x r- 1 / 2 w d ^ l ((n + 1 )r 2 ) , (7) 

where 7 = (d — 3)/2 is a dimension parameter for the Coulomb system that van- 
ishes in the three-dimensional case. The constant -£Q, n ,A is selected to preserve the 
normalization. The restrictions on this correspondence may be expressed as 

D = 2d-2-2\ , (8) 

N = 2n-2 + \ , (9) 

L = 21 + \ , (10) 

where it can be seen from flTU| ) that A, which gives an extra degree of freedom in the 
mapping, has to be integer valued. For d — 3, A may equal zero or one, yielding 
oscillator dimensions of D = 4, 2. This is Schrodinger's original result. Equation (|) 
shows there is no such correspondence between the physically interesting cases of the 
D = 3 oscillator system and the d = 3 Coulomb system. The oscillator is limited to 
even dimensions D only. 

To circumvent these dimensional restrictions on the correspondence between the 
exact systems, we broaden the class of systems considered to include ones with ana- 



lytical modifications of the type introduced above for interactions between particles in 
a trap. So, whereas there is no natural mapping between the single-particle Coulomb 
and oscillator systems, a map may exist between a trap with several fermions and the 
exact Coulomb system. 

We allow for an integral shift J in the oscillator dimension via the definition and 
requirement D* = D+J > 1. If we also define a dimension parameter V* = (£)*— 3)/2, 
which vanishes for D* = 3, then with the choice of effective potential 

h 2 (L* + T*)(L* + T* + 1) - (L + T)(L + T + 1) 



V KVV (r) 



1m r 2 

+ hu (N-N*+T-r*) (11) 



we obtain a differential equation with analytical solutions Wrj*,N*.L*( r )- The values 
of the parameters A, /, and J are restricted if normalizability and orthogonality 



are desired ||18|| . With this broader class of oscillator radial systems, the case of a 
mapping from the D* = 3 oscillator to the d = 3 Coulomb system becomes possible. 
It is obtained by setting 

A- I = X-l/2 , (12) 

and since A can only take values zero or one, the mapping requires a nonzero defect 
A in the oscillator. Explicitly, the relationship and its constraints are 

WV*,L*(r) = iW/2r- V Vw(^ 2 ) , (13) 

iV* = 2n-3/2 , (14) 

L* = 2/ + 1/2 . (15) 



Condition (13) ensures that the entire stack of states of the one system maps across 
to the entire stack for the other, with the lowest states in the stacks identified with 
each other, the second lowest states with each other, and so on. 

This mapping is not the only one possible between the three-dimensional sys- 
tems. An alternative method involves allowing for an analytical modification in the 
Coulomb system instead of in the oscillator and then following methods similar to 
those leading to flT3|). Such an analytical quantum defect may be introduced by shift- 
ing the principal quantum number to give eigenenergies according to the well-known 

6 



Rydberg formula E Q /(n*) 2 ||24)| , where E is the ground state energy of the Coulomb 
system and n* = n — 5, with 5 = 5(n, I). Analytical solutions exist for this system, and 
their applications include the study of highly excited valence electrons [|, |, |, §, [7|. 
Since the Rydberg formula models the spectra of multi-electron atoms, this option 
establishes a mapping between a single-particle trap and a multiparticle atom such as 
an alkali-metal atom. A more general third option involves analytical defects in both 
systems and provides a mapping from a multiparticle trap to a multi-electron atom. 
This most general form of the mapping incorporates in special cases the bosonic sec- 
tors, the fermionic sectors, and the quantum-defect sectors for various dimensions. 
For example, in a mapping between the three-dimensional cases, the generalization 
of ( |12"D can take the form 

A-I = 2(5-i) + X-l/2 , (16) 

where i is the analogue for the Coulomb problem of / for the oscillator. 

Unlike the defects A for traps, the defects S have been measured for many atoms 
||20|| . They depend on I but are asymptotically independent of n for large n. This 
feature is attractive as it allows approximate orthogonality of the solutions. It would 
be equally attractive if it could be established for the oscillator system. Indeed, if the 
oscillator defects A could be measured, it would be amusing to know if the constraint 
(P3|) on the radial mappings between the three-dimensional systems is consistent with 
the known atomic defects 5. 
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